Electron spin quantum beats in positively charged quantum dots: nuclear field effects 
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We have studied the electron spin coherence in an ensemble of positively charged InAs/GaAs 
quantum dots. In a transverse magnetic field, we show that two main contributions must be taken 
into account to explain the damping of the circular polarization oscillations. The first one is due to 
the nuclear field fluctuations from dot to dot experienced by the electron spin. The second one is 
due to the dispersion of the transverse electron Lande g-factor, due to the inherent inhomogeneity 
of the system, and leads to a field dependent contribution to the damping. We have developed a 
model taking into account both contributions, which is in good agreement with the experimental 
data. This enables us to extract the pure contribution to dephasing due to the nuclei. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

A single carrier spin in a single quantum dot repre- 
sents a potential candidate for q-bit implementation in 
a solid state environment, in view of applications in the 
fields of spintronics and quantum information processing 
[E H H H- Quantum dots (QDs) are indeed attractive 
with respect to the general criteria for quantum com- 
puters [5j], since long spin relaxation times T\ have been 
measured in neutral QDs, in the millisecond range for 
longitudinal magnetic field between 4 and 8 T pi, 0]. 
This is made possible due to the inhibition of the ran- 
dom interactions leading to spin relaxation and decoher- 
ence in bulk Q and quantum well semiconductor struc- 
tures. However, one of the major difficulties to achieve 
quantum manipulations is decoherence due to interac- 
tions with an uncontrollable environment. Long spin de- 
coherence times Ti are indeed demanded in order to be 
able to achieve enough quantum manipulations during 
the spin lifetime Ti, and this requirement turns out to 
be the most stringent one. In principle, the dephasing 
rate due to spin-orbit coupling may reach values of the 
order of T\ [9(. Recent experiments have demonstrated 
a trend in this direction, where Ti of the order of 1 fis 
could be measured at low temperature, allowing possible 
applications based on robust quantum coherence within 
an ensemble of dots [l(| E| ■ 

In quantum dots, two dephasing processes are still ef- 
ficient, and even enhanced by confinement: the exchange 
interaction between charges, and the hyperfine interac- 
tion. The first one is responsible for various phenom- 
ena, such as e. g. optically active exciton state split- 
ting, neutral exciton quantum beats and the appearance 
of negative circular polarization for negatively charged 
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excitons under non resonant excitation pjl, [IU HH, [H[ • 
The second one, the hyperfine interaction of localized 
electrons with the QD nuclei, leads to very efficient spin 
dephasing [l6|, [Tt], IU, This interaction effect has 

been observed previously in p-doped InAs/GaAs QDs 
for the ground state of the positively charged exciton 
(also called trion) X + by time resolved photolumines- 
cence (TRPL)[19|], and more recently in n-doped QDs 
for the resident electron using time resolved Kerr pump- 
probe spectroscopy [ll|, |2(j. Three distinct time scales 
are relevant in describing the electron-nuclei spin system 
evolution [l6l |: the first time corresponds to the electron 
spin precession around the frozen nuclear field fluctua- 
tions due to the QD nuclei (the typical dephasing time is 
of the order of 1 ns for GaAs QDs containing 10 5 nuclei); 
the second one is controlled by nuclear spin precession 
in the hyperfine field of the localized electron (the typ- 
ical time is of the order of 1 fis), and the last one is 
the nuclear spin relaxation due to dipole-dipole interac- 
tion with nuclei in the vicinity of the QDs (the typical 
time is of the order of 100 /is). During the first two 
stages, the coherence of the electron-nuclear spin system 
is preserved, while during the last one it is not, since the 
dipole-dipole interaction docs not conserve the total nu- 
clear spin. In ref.[l(| and it is essentially the first 
dephasing stage which is observed, leading to an estimate 
of the electron spin ensemble dephasing time T 2 * in the 
inhomogeneous nuclear field. The longer spin coherence 
times measured in ref. [ljj H3 are based respectively 
on spin echo technique or mode locking of electron spin 
coherences to suppress the hyperfine induced dephasing, 
and are presumably limited by the nuclear spin dephasing 
time. 

In this work, we investigate the first stage of electron 
spin dephasing in an ensemble of InAs/GaAs QDs in 
the presence of an in plane magnetic field. To achieve 
this aim, we study the electron spin coherence (ESC) 
in time resolved photoluminescence spectroscopy of pos- 
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itively charged excitons X + . Exciting QDs containing a 
single doping hole with an optical pulse results in the for- 
mation of a trion X + which, in its ground state, consists 
of a hole spin singlet in the highest valence states and a 
single electron spin in the lowest conduction orbital state: 

h4)-^(|i-i>-|4!»H4) w 

where |±§), (±3) represent respectively the projection 
of the heavy hole and conduction electron angular mo- 
mentum on the quantification axis Oz, taken normal to 
the sample surface. Hence, the electron-hole exchange 
interaction which is efficient within neutral QD excitons 
[3, HH is cancelled out [H, [22j • The time and polariza- 
tion resolved photoluminescence (TRPL) signal is thus a 
direct probe of the unpaired electron spin dynamics dur- 
ing the radiative time. Note that in this approach, when 
exciting the dots within the wetting layer (WL), the holes 
lose thei r spin p olarization before they are captured by 
the QDs [la, |23[ due to efficient spin relaxation in the WL 
[HI]. In addition, the spin coherence of localized holes 
that might be generated by the laser pulse, in a symmet- 
rical way as for negatively charged dots [23, HH, [2f| , will 
be cancelled out after the formation of the hole singlet 
in the trion ground state. Finally, the hyperfine interac- 
tion of the resident hole spin with the nuclei is negligible, 
due to the p-symmetry of the periodic part of the Bloch 
function [26|, [27| . Note also that, for n-doped QDs, the 
Kerr resonant pump-probe approach is very appropriate 
to study long term resident electron spin evolution, par- 
ticularly after the radiative recombination 11, 20]; how- 
ever, the response at short time delay is more complex 
to analyze using this technique: both the dots with a 
photogenerated trion and the dots with a single carrier 
bearing an optically generated spin coherence contribute 
to the Kerr probe signal, leading to a complex interfer- 
ence pattern under transverse magnetic field (due to non 
zero transverse g- factor of the hole in QDs [2l|). The 
TRPL experiment is then well adapted to make an accu- 
rate description of the unpaired electron spin coherence 
of X + during the first step of its dephasing by nuclear 
spins. As we observe the average electron spin (S z (t)) in 
an ensemble of dots, the decay of the oscillation ampli- 
tude can be described by a characteristic spin dephasing 
time T 2 * ; in a classical view, this is the decoherence time 
of the spin ensemble during the precession around the ap- 
plied transverse magnetic field, taking into account the 
different inhomogeneities of the sample. 

We have developed a theoretical model to describe the 
quantum beats observed experimentally. We observe that 
two contributions participate to the spin dephasing: the 
first one is due to the nuclear field fluctuations. However, 
this contribution alone cannot explain the observed de- 
phasing: the second contribution is due to static Lande 
g- factor fluctuations from dot to dot [2(| l25| - Finally, 
variations of QD doping level have also to be taken into 
account. By fitting our theoretical model to the experi- 
mental data, we can extract the contribution to the de- 



phasing due to the hyperfine interaction alone. The hy- 
perfine interaction is also responsible for the dephasing in 
a single quantum dot, when averaging over a large num- 
ber of nuclear field fluctuations. This is the case when 
the signal integration time is much larger than nuclear 
spin relaxation time due to the dipole-dipole interaction. 



II. SAMPLES AND EXPERIMENTAL SET-UP 

The samples used for this work consist of 10 planes 
of self assembled InAs QDs grown by Stranski- Kras- 
tanow method, and separated by 30 nm of GaAs. A 
beryllium delta doping layer is located 15nm below each 
wetting layer. Several samples have been investigated 
with different nominal doping corresponding roughly to 
one hole per dot on average. We show here the results 
for two typical samples : sample 1 with a nominal dop- 
ing of Na = 5 x 10 10 cm~ 2 , and sample 2 with a nom- 
inal doping of Na = 15 x 10 10 cm -2 . The samples are 
mounted in a cryostat cooled at a temperature of 15K. 
The experiment is performed in Voigt geometry, the ex- 
ternal transverse magnetic field B = B x e x being oriented 
along the Ox axis in the QDs plane. The QDs are ex- 
cited with 1.5ps pulses from a mode- locked Ti-sapphirc 
laser with a repetition rate of 80 MHz. Circularly po- 
larized light (er+) propagating along the growth axis Oz 
generates a coherent superposition of \X + . +1/2 > x and 
\X + , —1/2 > x trion eigenstates. The excitation beam is 
focused on a spot size of 100/Ltm diameter with an av- 
erage power of lmW. We did not observe any change 
in the PL circular polarization dynamics for excitation 
powers up to 5mW. The PL intensities co-polarized (7 + ) 
and counter-polarized with the excitation laser are 
dispersed by a monochromator and then recorded us- 
ing a SI photocathode Hamamatsu Streak Camera with 
an overall time-resolution of 30 ps. We measure the 
circular polarization degree of the photoluminescence 
P c = (I + — + I~) , which corresponds directly 

to the electron average spin component (S z (t)) = —P c /2 
along Oz @]. 

III. RESULTS AND DISCUSSIONS 

We present first the circular polarization (P c ) dynam- 
ics when increasing the magnetic field B x up to 750mT. 
The excitation energy (E ex = lAAeV ) cor resp onds to 
the lowest states of the wetting layer (WL) [23]. Under 
these excitation conditions we assume, consistently with 
previous experiments, that the electron maintains its spin 
orientation during the capture in QDs [23|. For sample 
1, figure [1] presents the time evolution of the circular po- 
larization P c (t) for three different magnetic fields B x . At 
zero magnetic field, we find positive circular polarization 
(not shown) confirming that the QDs are p-doped [l9| . 
We recall that for neutral QDs the anisotropic exchange 
interaction yields linearly polarized neutral exciton eigen- 



3 



o 

'•*—> 

03 

N 
"i_ 
_03 

O 
Q_ 

30 
O 20 
O 10 





Sample 1 >3o| 

~10- 
LU 

< Ob 



A 



250 500 750 

B v (mT) 



A A x 



B =750mT 




B =330mT 



B =190mT 







500 1000 

time (ps) 



o 

03 
N 
'l_ 

O 
□_ 

i_ 
ttJ 

O 

b 



30 
20 
10 




Sample 2 >3o 

^ 20 

< o 



\ 








250 500 750 

B v (mT) 




B =720mT 



B =320mT 



B =140mT 



500 1000 

time (ps) 



FIG. 1: Sample 1 (N A = 5 x 10 10 cmT 2 ): time-resolved cir- 
cular polarization decay of the X + trion under a+ excitation 
for three different magnetic field strengths B x . Black dotted 
lines represent experimental curves, grey lines are theoretical 
curves (see text). Inset: Zeeman splitting as a function of B x 



states and under non resonant excitation, the circular po- 
larization value is lower than 3% [HI, H3| • For n-doped 
QDs, the eigenstates are circular, but the polarization 
rate is negative due to interplay of anisotropic exchange 
interaction and Pauli blocking [15, 23, 30, ^31(. Moreover, 
as we can observe luminescence under strictly resonant 
excitation (not shown here), we deduce that there are on 
average less than 2 holes per dot. We assume first that all 
the dots contain a single resident doping hole (this point 
will be discussed later). After relaxation, the optically 
excited and the doping holes will finally form a singlet in 
the trion ground state, which recombines radiatively. 

The oscillation period observed, corresponding to the 
X + ground state Zeeman splitting AE, decreases when 
the magnetic field increases. The inset in FigfT] shows 
the linear dependence: AE = g±ofisB x . We can 
thus find the average transverse electron Lande g-factor 
Iff -Lo I ~ 0.75 ± 0.05. Note that, as in quantum wells 
|32l . [H, Hi], [HI, the electron effective Lande g-factor is 
strongly anisotropic [36| . 



FIG. 2: Sample 2 (N A = 15 x 10 10 cm" 2 ): same as Fig[T]but 
increasing the doping level of the structure. 



From the curves P c (t) plotted in FigQ]we also see that 
despite the relatively high magnetic field value, the ESC 
decays with a typical time T 2 * of about 300 ps. Increas- 
ing the magnetic field also leads to an increase of the 
damping of the polarization oscillations. Furthermore, an 
unusual asymmetry of spin quantum beats is observed, 
the negative extrema of the circular polarization being 
smaller in absolute value than the positive ones, in con- 
trast with what is commonly observed [3, [13, H3, [H[. 
Comparing Fig[TJ and FigJ^l we see that this asymmetry 
increases with the doping level of the structure. 

To explain these observations, we discuss first the role 
of the interaction of the trion with nuclei. Secondly, we 
will explain the dependence of the beats damping on B x 
magnitude and the reason for the asymmetry. Finally, 
we will show the agreement between the experimental 
data and the different theoretical parameters used in this 
paper. 

As the correlation time of the nuclear field is much 
longer than the X + trion lifetime [T§], we start from 
the model of the " frozen" configuration of nuclear spins 
[l6j |. In a first step, we only consider the electron spin 
dephasing originating from the random distribution of 



4 



the nuclear hyperfine field in the dots. 

In order to determine the equation of motion of a spin 
in a dot, we adopt a smilar approch as in Ref [l6(. How- 
ever here, care has to be taken in order to take into ac- 
count the g-factor anisotropy. It is thus more convenient 
to use the precession vector picture rather than the ef- 
fective magnetic field picture. The Hamiltonian of an 
electron in a dot is: 

H = H N + H Z (2) 
H = v Q A * l*( R *)| 2 US + ( VBgB ext ).S (3) 

i 

where the first term is the usual hyperfine contact Hamil- 
tonian, and the second is the Zeeman hamiltonian. Here 
g is the electron g-factor tensor, Ij and S are writen in H 
units, and vq is the volume of the unit cell. The hyper- 
fine interaction constants Ai are all of the order of 50/ieV 
considering two atoms per elementary cell |39(. We can 
thus write: 

H = h£l.S = h(£l N + £l ext ).S (4) 

where £In = vq / h^2 i Ai ^ (Ili)\ 2 li is the precession 
vector operator and £l ex t — HBgBext/h- The quantum 
mechanical average of fl^v on a given nuclear configura- 
tion is: 




In the " frozen fluctuation" approach, Jljy evolves with 

time much slower than the electron average spin ^S(i)^. 

The electron sees a distribution of Jljv, whose magnitude 
and direction are randomly distributed over the QD en- 
semble, described by an isotropic gaussian probability 
density distribution function : 

where Ao N represents the dispersion of the precession 
vector. The fluctuation of $7 at can be readily obtained as- 
suming that the nuclear spin directions are independant 
from each others [HI, : 
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where Nl is the number of nuclei interacting with the 
electron in a dot and n is the number of atoms in the lat- 
tice unit cell. The symbol ((. . .)) stand for averaging on 
the different nuclear field configurations. It is thus possi- 
ble to define a characteristic dephasing time: Ta = Aq^, , 
(independently of the g-factor tensor). The equation of 

motion of an average spin /s(t)^ in a fixed magnetic field 

in a QD is then given by: 

- (S .n)n+ [S - (S .n)n}cos(fl t) 

+ [S - (S .n)n] x n sin(Q t) (9) 

where: ft = J7at + J^ezt, n = fl/ \Si\ is a unit vector, 
r2jv is the nuclear precession vector in a given QD as de- 
fined in (O and S is the initial average electron spin 
Qext = 9ef f^BB ex tlh represents the contribution of the 
external field to the total precession vector $7. In this ex- 
pression, an effective Lande g-factor has been introduced, 
which is defined by : <? 2 y^ — 9± sin 2 r/ + g?, cos 2 rj, where 
V = (e z , ~B e xt) is the angle of the external magnetic field 
with the Oz axis. The angle rf = (e Zl tt ext ) is given by 
rj = Arctan(^tan rj). 

Expression ([9]) is then averaged over the QD ensemble, 
taking into account the Gaussian variations of £In char- 
acterized by the equation §§§ 0j[ . The calculation is 
performed in the Appendix for an arbitrary orientation 
of the external magnetic field (the Ox axis is chosen in 
the plane defined by So and B ext ). We assume that no 
dynamical polarization of the nuclei occurs; so that the 
ensemble average of JIn is zero. For a pure transverse 
external field (Voigt configuration) this assumption is al- 
ways valid, since no dynamical polarization of the nuclei 
(i.e. Overhauser effect) can occur [44| . 

The average electron spin S(t) = ^/lS(i)\^ can be 

expressed as the sum of two contributions S(t) = S°° + 
Si(f), where S°° is time independant and Si(t) contains 
the oscillating contribution which damps as t increases. 
The expressions S°° and S\(t) are given respectively by 
the expressions (|A9l and ([AlO]) . 

For weak magnetic fields (B x < 200mT), taking Ta = 
500ps as deduced from the measurement at zero mag- 
netic field [l9] leads to satisfactory fits of the experimen- 
tal data. However, for larger magnetic fields, we will see 
below that the contribution of the nuclear field fluctu- 
ations is not sufficient to explain the observed damping 
(see FiglSJ). Moreover, the dephasing time Ta is field 
dependent, which contradicts the fact that, under trans- 
verse magnetic field, the nuclear field fluctuations are in 
principle insensitive to the applied magnetic field. 

In fact, the magnetic field dependent damping arises 
from variations of the electron g-factor over the QD en- 
semble, leading to a sprea ding of the Larmor frequen- 
cies with increasing B x 0, |25| | For instance, the origin 
of this inhomogeneity could come from the different dot 
sizes, or chemical repartition of In in the QDs. This 
contribution is the second source of spin dephasing for 
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FIG. 3: Theoretical curves from the model developped in the 
appendix showing the spin quantum beats in a weak trans- 
verse magnetic field (B x — 150mT)when the amplitude of 
the g-factor and the nuclear field fluctuations vary, (a) fixed 
nuclear field fluctuations (Ta = 500ps) and increasing g fluc- 
tuations, (b) fixed g-factor fluctuations (A g /g±o = 0.07) and 
increasing the nuclear field fluctuations (i.e. decreasing the 
parameter Ta- 



FIG. 4: Same as Fig[3]but for ; 
field (B x = 750mT) . 



stronger transverse magnetic 



tween the two cases. The fluctuations of the precession 
vector due to external magnetic field are now given by 
the gaussian distribution: 



an ensemble of quantum dots. Under the reasonable as- 
sumption that the nuclear field and g-factor fluctuations 
are uncorrelated, we have included the variations of the 
latter in the model of Merkulov et al. [l6| . For an arbi- 
trary magnetic field orientation, we did the simplifying 
assumption that the transverse and longitudinal fluctua- 
tions (5g± and 6g\\ respectively) are strongly correlated, 
according to 



5g\ 



9± 

9\\ 



gxo 
9\\o 



(10) 



where <?j_(||)o is the average <7_l(||) value. Under this as- 
sumption, only the amplitude of fl ex t fluctuates from dot 
to dot, but not its direction, which is given by the an- 
gle rj' = arctan(^- tan n). This approach yields exact 

results for pure transverse or pure longitudinal magnetic 
fields. For oblique fields, it gives some interpolation be- 
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sin 2 rj 



and Aq, defined by A 2 ^ = Ag^ sin 2 r/ + 
Ag 2 cos 2 n, characterizes the fluctuation of the precession 
vector length. Note this corresponds to a gaussian proba- 
bility distribution characterized by the parameter A 9e// , 
and the average value g e ff,o- 

Then expression (|9]) can be averaged taking both gaus- 
sian variations for Bn and g±. The new equations for 
S°° and Si(t) are now given respectively by the expres- 
sion (|A9j) taken for g e ff,o and (|A22j) . from which the time 
resolved circular polarization can be deduced. 

Let us now explain the different effects of the nuclear 
field fluctuations and the g-factor fluctuations. The ob- 
tained theoretical curves derived from the model devel- 
opped in the appendix are shown in FigJ3] for B x — 
150mT and in Figg]for B x = 750mT. In FigETa) and 
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FigHJa), the nuclear field fluctuations are kept constant 
(T\ = 500ps), while the g- factor fluctuations increase. 
In FigCJb) and FigHJb) A g /g 0± is kept constant at 0.07 
while the nuclear fluctuations increase. From the com- 
parison of both figures it is clearly seen that, at weak 
transverse magnetic field, the damping of the oscillations 
is mostly determined by the nuclear field fluctuations, 
while both contributions are necessary under larger mag- 
netic field using realistic sets of parameters. 

Note that the curve in FigJ3Ja) for g- factor fluctuations 
A g /g±o = allows us to extract the contribution of the 
hyperfine interaction (the equation of the spin quantum 
beats is given by the expressions (|A9[) plus (|A10|) ). We 
clearly see that the nuclear field fluctuations lead to a first 
contribution to damping of the P c oscillations even for the 
strongest B x applied. The transverse magnetic field B x is 
responsible for a large amplitude of spin quantum beats 
and yields a precise measurement of the spin decoherence 
time T 2 * . The characteristic dephasing time appearing in 
the model at zero magnetic field is Ta = A^ which 
typical value is here Ta = 500ps. The characteristic 
time for the decay of the oscillations amplitude is 2Ta 
(see equation (|A10I0 . 

The variations of P c (t) according to the model includ- 
ing the g- factor fluctuations are displayed on figure for 
A g /g±o = 0.07 and A g /g± = for comparison. We 
see that, as expected, the g factor fluctuations lead to 
an increase of the dephasing. The new characteristic de- 
phasing time T 2 for an ensemble of dots is now given 
by: 




where Ag = h/ g e f ffJ-sTA- The variations of T 2 /Ta 
with the magnetic field are displayed on figure [5] (bold 
line). When the magnetic field is much stronger than the 
nuclear field fluctuations (B » Ag -9eff,o/ A Se// ) , these 
variations can be approximated by : 

= t^^ Ta (13) 

The corresponding dotted line in figure which is the 
expression taken in ref. [20l |. is given for comparison. 

At this stage, the model correctly describes the period 
and the damping of the electron spin quantum beats. 
For transverse magnetic field comparable to the nuclear 
field fluctuation the beats are non symetrical with re- 
spect to P c = 0. However, in the investigated range, B x 
is stronger than As, which result in symetrical P c (t) os- 
cillations. Thus, the observed asymmetry remains to be 
explained. 

We believe that the origin of this asymmetry lies in 
the QDs charge fluctuations. Besides dots doped with 
a single resident hole, some are neutral and some con- 
tain two resident holes. So that neutral excitons (X°) 
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FIG. 5: Sample 1: comparison of experimental (dotted line) 
and theoritical curves of time resolved circular polarization 
for B x = 750mT. The theoretical curves (grey line) is given 
with (A g /go± = 0.07) or without (A g /go± — 0) g-factor fluc- 
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FIG. 6: Dephasing time as function of the relative magnetic 
field. Dotted line shows the exponential model used in ref 
[20I ]. Bold line represent our calculations (see equation I12H 

and doubly charged excitons (X 2+ ) are detected. The 
anisotropic exchange interaction (AEI) between electron 
and hole in X° and X 2+ would lead in principle to beats 
of P c (t) for these QDs even for B x = [lj,ll5|. However, 
the latter quickly damp due to the strong disp ersion of 
the exchange energy from dot to dot [U, |U |42|, |43| . 
For X° and X 2+ , the spin polarization decay can be 
described by: S(i) = Soe~ t / T °> where T* is an ensem- 
ble inhomogeneous dephasing time. This expression is 
valid provided that the dispersion of the AEI energy is of 
same order of magnitude as its average value. Defining 
a as the relative number of dots with charge state other 
than one, the time evolution of the circular polarization 
isgiven by: P c (t) = -2S ((1 - a)S(t)/S + oe"*/^) 
[451 ]. The theoretical curves of P c (t) on Fig[5](grey lines) 
are displayed using the parameters Ta = 500ps, a = 0.4 
and T* = AOOps, with (A g / 5o _L = 0.07) and without 
(Ag/gox — 0) g-factor fluctuations and are compared to 
the experimental curve (dotted line). 

To summarize, the only adjustable parameters of the 
model are A g /g± , T*, and the weight a of QDs with 
or 2 holes (depending on the doping distribution), all 
independent of the magnetic field. The nuclear field flue- 



7 



tuations in the QDs ensemble are determined from a fit 
of the circular polarization without magnetic field p^ |. 
and the average transverse g-factor g±o from the period 
of the oscillations (see fig.l). Note that A g /g±o = 0.07 
is consistent with previous studies (36j . 

The curves on Figures [Hand [2] are fitted using a unique 
set of parameters A g /g± , T*, changing only the mag- 
netic field. The parameter a takes into account that 
sample 2 (Na — 15 x 10 10 cm -2 ) is doped more than 
sample 1 (N A = 5 x 10 10 cm~ 2 ). We take a = 0.4 for 
sample 1 and 0.6 for sample 2. We observe a good agree- 
ment between theory and experiment for several applied 
transverse magnetic fields. 

In conclusion, we have developed a model to describe 
electron spin dynamics in an ensemble of p-doped QDs 
under transverse magnetic fields. We have shown that: 
(i) the hyperfine interaction, (ii) g-factor fluctuations 
from dot to dot, and (iii) QD charge variations, have 
to be taken into account to explain the time evolution 
of the emited circular polarization P c . The experiments 
presented here demonstrate the key role played by the 
nuclear field fluctuations on the electron spin coherence 
dynamics in QDs and we found a dephasing time of about 
1 ns for this contribution. Note that in single QD exper- 
iments, only the contribution due to the interaction with 
the nuclei should be observed. 
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APPENDIX A: THEORETICAL MODEL 



We use in this appendix ft 



Ac, ' Z 



as reduced parameters in units of the nuclear contribu- 
tion to the precession vector fluctuation. The parameter 
z is taken algebric and we recall that ft = tt ex t + £In — 
fl.n. 

Note that the relation: 



COO 000 

2^ + 2^ = 1 

Sox S z 



(A3) 



is fulfilled (see ref. [Ty|) whatever the value of ft is. 

The second contribution contains the oscillating contri- 
butions which damp with a typical time 2Ta- We define 
the following integrals: 



„t 2 /4-,3 2 



2 e 2/3z - 1 
e z — — cos (zt) dz (A4) 



2ftz 



/(/3,r) = — (co8(Pt)I 2 (J3,t) 
+ T -sin{PT)h{p.T) - e-? 2 7 2 (0,t)) 
using the the following auxilliary integrals: 



V7T 



-oo Z 2 +T 2 /A 



dz 



-dz 



(A5) 



(A6) 



(A7) 



, Z 2 + T 2 /A 

We can thus rewrite the different components of S 



as: 



2f3 2 



(l-/(/3,0)) 



We first take into account only the fluctuations of the 
hyperfine field represented by the equation |6]). 

The calculation is performed in such a (OXYZ) frame 
that the precession vector fl ex t is along the OZ axis (r/ = 
(e z ,ez)) and the initial electron spin axis S(0) is in the 
(XOZ) plane, i.e. OX belongs to the (S(0), Q e xt) plane. 

Averaging the equation §§§ of the spin motion we can 
find the evolution of the mean spin as a function of time, 
we obtain: 



C X 

and thus: 

ooo _ Sqx 

Sx ~ w 



= 



(A8) 



Soz 



(/3 2 -1 +!(/?, 0)) 



L(l 2 (/3,0)-e-^/ 2 (Q,0)) 



2ft 



S(i) = S°° + S 1 (i) 



(Al) 



The first contribution is time independant. It repre- 
sents the asymptotic limit for S(t) when t — > oo. The 
different values of S(t) in the (OXYZ) frame are: 



qoo ^OX 

bx - w 



1 - 



1 



-(* 2 +/3 2 ) 



2-, e 



2/3z _ 1 



= 



2ftz 



-dz 



(A2) 



Sf = 



1 



(A9) 



ft 2 -l + ^(h(P,0)-e-P h(0,0) 



The oscillating contributions Si(t) are given by the 
following expressions: 



/3 2 



2ftz 



dz 



8 



Six® = S 0X e^ 2 /4 |(i _ ^)cos(0t) ~ ^sin(pr) + ^/(Ar)} 



Siy(t) - -W^ 4 (1 



^)sin(0T) 



2(3 



cos(/3t) 



(A10) 



S lz (t) = ^e- 2 /4 



{co S (/3r)-/(/3,r)} 



Note that the following relation holds: 



d_ 

dt 



2S lx (t) S lz (t) 



5, 



OA- 



S', 



oz 



= 2o (Srt ^rW (A11) 



OA 



For the integral 1(0, 0) the series expansion can be given: 



WO) 



E 



p 2 p 



(A12) 



From this expression we can calculate 1(0, 0) and find the 
same limit as in Ref.pjj]: 



S 



A 



OA, 



tS* v' — Oj s 



(A13) 



For small magnetic field (0 « 1), the general expres- 
sion of S(t) is, up to second order with respect to (3 : 

Sx(t) = ?f jl + 2e - 2 / 4 (l-^)}+0(/? 2 ) 



s Y (t) = o 



(A14) 



Sz(t)= S -f {l + 2e-^(l-^)} 



0((3 2 ) 



An other limit case is found at high magnetic field 
(0 » 1). The expression for S°° is: 



S 



A 



0; S£?«0 



oz 



and the expression for Si(t) is: 

Six(t) = S QX e- T2/4 cos((3T) 

S 1Y (t) = -S QY e- r2/4 sin([3T) 

Siz(t) = 
and finally the expression for S(t): 

5x(t) =S QX e- T2 / A cos(f3T) 



(A15) 



Sy(i) = -S , oy e~ r2 / 4 sin(^r) 



(A16) 



S z (t) = 5, 



oz 



Note that, at long time delay, we have Si(i) — > and 
S(t) -> 5*°° when r -► oo. 

Let now discuss the averaging of equations (|A9[) and 
(|A10() if we take into account the g- factor fluctuations. 
We introduce the effective g-factor by the expression: 



ex I 

h 



where: 



2 2-2,22 

9 e ff = 9±sm rj + g^cos r) 



(A17) 



(A18) 



and r\ = (e z ,B ext ). 

The g-factor fluctuations lead to fl ex t variations given 
by the expression (jTTJ) where: 



An = 

n 

Using the expression j3 = fl ex t/A.Q N we obtain: 
We assume in our model that: 



(3 



__9ef£ 
9eff 



« 1 



(A19) 



(A20) 



(A21) 



Moreover, if > it (i.e. B ext » h/2g e fffj,BT&), the 
envelope of the average spin oscillations decays much 
slower than the period of these oscillations. It follows 
that, in order to average S(t) f (|A9[) and (|A10[1 ) over the 
g-factor fluctuations, we can use the following approxi- 
mation: S°° » S°°)p = i3 . For Si(t) we obtain the ap- 
proximate expression: 



9 



Six(t) = Soxe-^+^f: 



1 

Wo 



J 2 (A),t) 

2/?o 



cos([3 t) 



2/3o 



sin(f3 T) \ - S ox e " 4 



e- Ql 2 (O,T) 
4/3 3 



2A 2 



o / 2 Po 



(A22) 



Po 



1 - 



2/3o 



cos(/3 t) 



The calculation frame (OXFZ) is obtained from the 
initial one (Oxyz) by performing two rotations. One with 
Oy axis with an angle n' : Roy(f]')(Oxyz) = (Ox'y'z 1 ) 
such that Oz' / /n ext . 

The second with Oz' axis and angle <j>: 
R 0y ((j))(Ox'y'z') = (OXYZ) such that the plane 
(XOZ) contains So- We can thus obtain the general 
expression of the electron spin evolution in the frame 
(Oxyz) by: 



S v (t) 



COST) COSIp 

sincf) 
—sinrj 1 cose, 



—cost] smq 

COS<f) 

sinrj' sincf) 



smvj 


COST]' 



fSx(t)\ 

Sy(t) 

\Sz(t) 



r h(J3 ,T) 
2/3 2 



sin(Por) > + Soze 



zle-g° a J 2 (0,T) 
2/?3 



For a pure transverse magnetic field n = n' = tt/2 and 
= so the transformation takes the simple form: 




Sz(t) 

s Y (t) 

Sx(t), 



(A23) 



In addition, as here So = So z e z, we have: S x (t) = 
S y (t) = and S z (t) = —Sx(t). Some example are given 
in figure 1, figure 4 and figure 5. 
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